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Abstract 

We study the longitudinal excitations of quantum antiferromagnets on a 
triangular lattice by a recently proposed microscopic many-body approach 
based on magnon-density waves. We calculate the full longitudinal excita- 
tion spectra of the antiferromagnetic Heisenberg model for a general spin 
quantum number in the isotropic limit. Similar to the square lattice model, 
we find that, at the center of the first hexagonal Brillouin zone T(q = 0) 
and at the magnetic ordering wavevectors ±[Q = (47r/3,0)], the excita- 
tion spectra become gapless in the thermodynamic limit, due to the slow, 
logarithmic divergence of the structure factor. However, these longitudi- 
nal modes on two-dimensional models may be considered as quasi-gapped, 
as any finite- size effect or small anisotropy will induce a large energy gap, 
when compared with the counterpart of the transverse spin-wave excitations. 
We also discuss a possible second longitudinal mode in the triangular lattice 
model due to the noncollinear nature of its magnetic order. 

PACS numbers: 75.10Jm, 75.30.DS, 75.50.Ee 
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1 Introduction 



Anderson's spin-wave theory (SWT) provides a good description for the low tem- 
perature properties of many two-dimensional (2D) and three-dimensional (3D) 
quantum antiferromagnetic systems on a bipartite lattice [1], including quantum 
corrections to the classical Neel ground state with two alternating sublattices and 
the doublet transverse spin-wave excitations. The quantum antiferromagnetic sys- 
tems on a triangular lattice, however, are the prototypal system with the frustrated 
spin alignments in the classical picture and have been under intensive study over 
the last few decades for fundamentally different type of ground states [2-4] . It is 
now firmly established by various methods [5, 6], including a SWT based one 
three- sublattices [7-17], that the ground state of the antierromagnetic Heisen- 
berg model has the long-ranged noncollinear order of the 120° magnetic three- 
sublattice structure with three transverse, gapless spin-wave excitations. 

Most isotropic antiferromagnets in one-dimension (ID) with low quantum 
spin numbers do not show Neel-like long-ranged order in the ground state due to 
the strong quantum fluctuations. The low-lying excitation states are also different 
from the 2D and 3D counterparts. In particular, by the exact solutions using the 
Bethe ansatz [18], the low-lying excitation states of the ID spin-1/2 Heisenberg 
model have been shown corresponding to the spin-1/2 object (spinons) where the 
spin-wave-like excited states are the triplet states of spinons which always appear 
in pairs [19], contrast to the doublet excitation states by SWT; the excitation states 
of spin-1 Heisenberg model for the linear chain, including the longitudinal one, 
have an excitation gap above the singlet ground state, first predicted by Haldane 
[20]. These theoretical predictions have later been confirmed in the antiferromag- 
netic compound KCuF 3 for spin-1/2 chains [21] and CsNiCl 3 for spin-1 chains 
[22] by neutron- scattering experiments. One interesting remaining question is 
whether or not there exist longitudinal excitations in quantum antiferromagnetic 
systems with long-ranged classical order in low temperature, as such modes will 
represent the oscillations in the magnitude of the long-ranged order parameter. 
The answer is affirmative. There is now ample evidence of the longitudinal ex- 
citation states in various quasi- Id structures with the Neel-like long-ranged order 
at low temperature, including the hexagonal ABX 3 -type antiferromagnets with 
both spin quantum number s = 1 (CsNiCl 3 and RbNiCl 3 ) [23, 24] and s = 5/2 
(CsMnI 3 ) [25, 26] and the tetragonal structure of KCuF 3 with s = 1/2 [27]. More 
recently, a longitudinal mode was also observed in the dimerized antiferromag- 
netic compound TlCuCl 3 under pressure with a long -ranged Neel order [28]. To 
our knowledge, no observation of longitudinal modes in any 2D or quasi-2d an- 
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tiferromagnets has been reported yet. Clearly, such longitudinal modes, which 
correspond to the oscillations in the magnitude of the magnetic order parame- 
ter, are beyond the usual SWT which predicts only the transverse spin-wave ex- 
citations, usually referred to as quasiparticle magnons in the antiferromagnetic 
systems. There are several theoretical investigations in these longitudinal modes 
using the field theory approach, such as a simplified version of Haldane's the- 
ory for the spin-1 systems [29, 30] or the sine-Gordon theory for the spin- 1/2 
systems [31, 32], and both treating the inter-chain couplings as perturbations. 
A phenomenological field theory approach focusing on the spin frustrations of 
the hexagonal lattice of the ABX 3 -type antiferromagnetic systems has also been 
made [33]. We recently proposed a microscopic many-body theory based on the 
magnon-density waves for the longitudinal excitations of spin-s quantum antifer- 
romagnetic systems, using the original bipartite spin-lattice Hamiltonians [34, 35]. 
The basic physics in our analysis follows Feynmann's theory on the low-lying ex- 
cited states of the helium-4 superfluid [36, 37]: the longitudinal excitation states 
in a quantum antiferromagnet with a Neel-like order are identified as the collec- 
tive modes of the magnon-density waves, which represent the fluctuations in the 
long-range order and are supported by the interactions between magnons; these 
longitudinal excitation states are constructed by the magnon density operator s z 
in contrast to the transverse spin-flip operator s ± of the magnon states [38]. Our 
numerical results [38] for the energy gap values at the magnetic wave vector are in 
good agreement with the experiments for the energy gap observed in the tetrago- 
nal structure of KCuF 3 with s — 1/2 [27] . We hope that more experimental results 
for the energy spectra at other wavevectors will be available for comparison. 

In this article, we extend our microscopic approach to study the longitudinal 
modes in quantum antiferromagnets on a triangular lattice where the magnetic or- 
der is noncollinear hence there are possible more than one longitudinal modes. We 
employ the approximate SWT ground state in our calculations and find that one of 
the longitudinal modes is gapless for the isotropic models in the thermodynamic 
limit, due to the slow divergence of the structure factor, but any finite size effect 
or small anisotropy will induce a large energy gap at the magnetic wavevectors 
when compared with the counterpart of the transverse spin-wave spectra, similar 
to the collinear square lattice model. We organize this article as follows. For com- 
pleteness, we briefly outline the main results of spin- wave theory for the triangular 
lattice model in Sec. 2, using the one-boson approach. We then apply our micro- 
scopic theory for the longitudinal excitations in Sec. 3, using the approximated 
ground state from SWT. We find that a large energy gap can be induced by a very 
tiny anisotropy, similar to the square lattice model. We also discuss the possible 



3 



other longitudinal mode of the triangular lattice and the possible extension of our 
calculations to the more realistic models in the quasi- Id systems such as those 
hexagonal ABX 3 -type antiferromagnets mentioned above in the last section. 

2 Spin-wave theory for triangular lattice models 

The classical ground state of the antiferromagnetic Heisenberg model on a trian- 
gular lattice consists of three alternating sublattices with spins on each sublattice 
align at an angle of 120° to the other two sublattices. The distance between the two 
nearest-neighbor spins of the same sublattice is \/3 of the lattice spacing which is 
taken as unity in this article. The spin-s Heisenberg Hamiltonian is given by 

H = jJ2^-S„ (1) 

where J(> 0) is the coupling parameter and the sum on runs over all the 
nearest-neighbor pairs of the triangular lattice once. Following Singh and Huse 
[9] and Miyake [10] , it is convenient to transform the Hamiltonian of Eq. (1) by 
rotating the quantum projection axis of the spins along the classical direction in 
the xz-plme at the six z-sublattices surrounding the one j-sublattice. This trans- 
formation leads to the following rotated spin operators, 

Sf^SfcosCQ.rO+S'sinCQ.rO, 

S? 5?, (2) 

S'f->5 , fcos(Q.r i )-5 , fsiii(Q.r i ), 

for all i sites, where Q = (4.7r/3, 0) is the magnetic ordering wavevector at the 
corner of the hexagonal Brillouin zone of the triangular lattice (see Fig. la). After 
the rotation of Eq. (2), the spin-wave theory can be formulated in terms of only 
one set of bosons, rather than three sets originally employed [8]. The Hamiltonian 
operator of Eq. (1) after this transformation is given by 

H = JJ2 [ cos tQ • ( r * " r M S ! S j + SfSj) + ASfSf 

+ sm[Q-(T i -T j )](S?S?-S?S?)], (3) 

where we have also introduced an anisotropy parameter A(< 1) along the y-axis. 
We will see later that even a very tiny anisotropy, such as A = 1 - 1.5 x 10~ 4 , 
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will induce a large energy gap for the longitudinal excitation spectrum at the mag- 
netic wavevectors ±Q when compared with the counterpart of the spin- wave spec- 
tra. Using the conventional Holstein-Primakoff transformations, the Hamiltonian 
of Eq. (3) can be expressed in terms of boson operators and a] in a series 
in power of s. The classical result for the ground-state energy is given by the 
0(s 2 ) term. The linear terms (i.e., linear in and a\) cancel each other out. The 
quadratic terms are in O (s) and are retained in the linear SWT. The cubic terms 
(in 0(y/s)) and the quartic terms (in O(s )) have been treated as perturbations for 
the higher-order corrections to the linear SWT [10, 14, 39]. Therefore, the Hamil- 
tonian of the linear SWT is given by, keeping only the 0(s 2 ) and 0(s) terms and 
after a Fourier transofrmation for the boson operators with the Fourier component 
operators a q and al, 



H' = --JNs 2 + 3JsJ2 



A q a q a_ q — —Bq^a^a^q + a q a_ q ) 



(4) 



where A q and B q are defined by 



4z = 1 + ( A - '- ) T „ B q 



respectively and ^ q is defined as usual by 

^ = ^E eiq ' rp = ^( cos ^+ 2 



Qx v3 
cos — cos — q y , . 



(5) 



(6) 



with the summation over the nearest-neighbor index p and the coordination num- 
ber z = 6 for the triangular lattice. The Hamiltonian H' is diagonalized by the 
canonical Bogoliubov transformation, a q = u q a q + v q a_ q , 



H' = ~JNs(s + 1) + Y,S q (a\a q + h, 



(7) 



where, £ q = zJsco q /2, is the spin- wave excitation spectrum with the dimension- 
less spectrum co q given by 



0J n 



7g )(l + 2A 7g 



(8) 



We plot the first Brillouin zone (BZ) of the triangular lattice and the dimension- 
less energy spectrum co q of Eq. (8) in Fig. 1. This spin-wave spectrum has three 
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zero modes, one at the center of the zone T(q = 0), and the other two at the 
corners of the BZ, the ±Q ordering wavevectors. The magnitude of the spin- 
wave velocity near these points are different, with values of = 3 and 
v ±q — 3Js\/3/2\/2 respectively. 
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Figure 1 : (a) The hexagonal first Brillouin zone of a triangular lattice in reciprocal space. 
The coordinates of the labeled points are, T = (0,0), P = (27r/3,0), L = (tt,0), 
Q = (4vr/3,0), M = (7r, 7T/V3), K = (2vr/3,2vr/Vk) and O = (0,tt/\/3). (b) The 
intensity of the linear spin-wave energy spectrum (dimensionless) uj q of Eq. (8) for a 
triangular-lattice Heisenberg model at A = 1. 



The sublattice magnetization m is given by m = s — rj with the magnon den- 
sity r] defined as the ground-state expectation value of boson number operator 
1] = (a\ai)g with value rj = 0.261 for the isotropic model in the linear SWT. 
This magnon density of i] = 0.261 per lattice site is larger than the corresponding 
number of 0.197 for the square lattice model, indicating the larger quantum correc- 
tion in the triangular lattice model due to the frustration. Furthermore, we notice 
that the three sublattice magnetization are noncollinear with 120° between their 
directions, hence we may have two longitudinal modes in the triangular lattice 
model, one corresponding to the 2-component magnon-density fluctuations in our 
rotated spin basis discussed above, the other corresponding to the perpendicular 
x-component magnon-density fluctuations. In the following section, we will dis- 
cuss the longitudinal excitations in the ^-component magnon-density fluctuations, 
similar to the collinear square lattice case where there is only one longitudinal 
mode and leave the discussion of the other possible longitudinal mode to the last 
section. 
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3 The longitudinal mode of the ^-component magnon- 
density fluctuations in triangular lattice model 

As mentioned earlier, the longitudinal excitations in an antiferromagnetic system 
with a Neel-lrke long-ranged order correspond to the fluctuations in the order pa- 
rameter. We have proposed a microscopic approach based on the magnon-density 
waves and have applied to bipartite antiferromagnetic systems where the mag- 
netic order is collinear [38]. Here we extend our analysis to triangular lattice 
models using the one-boson approach for the approximation of the ground state 
and investigate the longitudinal mode of the ^-component magnon-density fluctu- 
ations. Briefly, the excitation state of the ^-component magnon-density waves is 
constructed by applying the magnon-density operator X q to the ground state \^f g ), 

\* e )=X q \* g ), (9) 

where the density operator X g , instead of the spin-raising and lowering operators 
s^ 1 for the transverse spin-wave excitations (i.e. magnons), is given by s z operator 
as 



l 

where the index / runs over all lattice sites and the condition q > ensures that 
the excited state |\l/ e ) is orthogonal to the ground state. We notice that in our 
rotated spin basis of Eq. (2), the three sublattice magnetizations all point in the 
same ^-direction. The longitudinal excitation spectrum is then given by 

where N(q) is given by the ground-state expectation value of a double commutator 

as 

N{q) = ±([X- q ,[H,X q ]]) g , (12) 
and the state normalization integral S(q) is the structure factor of the lattice model 

l,V 

In Eqs. (12) and (13), the notation (. . . ) g indicates the ground-state expectation. 
It is straightforward to derive the double commutator for triangular lattice by em- 
ploying the one- sublattice Hamiltonian of Eq. (3) after the rotation to obtain the 
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following general expression, 

N ® = w J ^ + 2A)(1 + riW 15 *}' + V 1 - 2A ^ - ^W s r +P h 



hp 



- sin[Q ■ (r, - r l+p )] [(S?S? +p ) g - (Sf S? +p > J 



(14) 



Now we apply the approximation of the linear SWT for the ground state in the 
expectations (. . . ) g and obtain 



N(q) = \s.jJ2 [(1 + 2A)(1 + lq )~g p + (1 - 2A)(1 - lq )~g' p 



(15) 



p 

7.1 



where the transverse correlation functions g r and g r are defined as 



both independence of index / due to the lattice translational symmetry. Their 
Fourier transformations are obtained as, using the approximation of the linear 
SWT for the ground state in the one-boson approach discussed in Sec. 2, 

g q =- , Bq =, g'=-(^Jb=-l), (17) 

and where A q and B q are as given by Eqs. (5). We obtain the numerical results 
at the isotropic point A = 1 as g p = 0.258 and g' = 0.034 for all the six nearest 
neighbors. As can be seen, N(q) is dominated by g p . 

Within the same approximation of the linear SWT for the ground state, the 
structure factor S(q) of Eq. (13) is obtained as 

<%) =V + Jf^2Vq'Vq+q' + Yl 9q'9q+q' , 08) 
q> q> 

where g q is as given by Eq. (17), rj is the magnon density as discussed before and 

"•^^Fir 1 '- <19) 

We notice that the structure factor of Eq. (18) involves the double magnon spec- 
trum function tu q cu q+q i in the integrals, indicating some effects of the interactions 
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between magnons has been included. We then obtain the longitudinal excitation 
spectrum E(q) of Eq. (11) using the approximations of Eqs. (15) and (18). By 
numerical evaluation, we notice that this spectrum of the longitudinal mode is the 
gapless in the thermodynamic limit, approaching zero E(q) — > at both q — > 
and q — > ±Q. Detailed numerical analysis shows that the gapless spectrum is 
due to the slow, logarithmic divergence in both the second and third terms in the 
structure factor S(q) of Eq. (18). More specifically, near T and ±Q, we find that 
S(q) oc — lng, and thus the excitation spectrum E(q) oc — 1/lng as q — > or 
q — > ±Q, but with different coefficients for the T point from that of the ±Q 
points. This later feature is reminiscent of the different spin-wave velocities at 
these points as discussed in the previous section. The major difference between 
the longitudinal and transverse modes near T and ±Q points will become apparent 
once we introduce the anisotropy as discussed in the followings. 

The logarithmic behaviors of the structure factor and of the energy spectrum 
of triangular lattice model is similar to that of the square lattice model studied 
earlier [34, 35, 38], where we state that the spectrum is quasi-gapped as any finite 
finite size effect or anisotropy will induce a large energy gap when compared with 
the counterparts of spin- wave spectrum. Here we consider the effect of anisotropy 
in the y-axis as given by Eq. (3). For a small value of A = 1 — 1.5 x 10~ 4 , we 
obtain energy gap value of 0.203zs J at the ordering wavevector ±Q, compared 
with the much smaller gap value of 0.0075zsJ of the corresponding spin-wave 
spectrum. More specifically, we find that the longitudinal energy gap value is 
proportional to l/[— m(l — A)], in contrast to the spin-wave gap which is pro- 
portional to y/1 — A, when A — > 1. We notice that the spectrum is still gapless 
at T due to the fact that the anisotropy A is introduced in the y-component in the 
Hamiltonian of Eq. (3) while the spins are ordered along the rotated ^-axis. In or- 
der to make further comparison between the longitudinal mode and the transverse 
spin- waves mode, we plot both the spectra with A = 1 — 1.5 x 1CT 4 in Fig. 2 
along the path (LMTKPQMO) of the BZ. The different gap values for the lon- 
gitudinal and transverse mode at both K and Q points can be clearly seen. Both 
the longitudinal and spin- wave spectra at T(q = 0) are still gapless where 7 q = 1 
for any value of A(< 1). Furthermore, there are slight different peak values at 
for the longitudinal spectrum but for the spin-wave spectrum the peak values are 
nearly equal. This is again similar to the case of the square lattice model dis- 
cussed earlier [38], indicating some effects of the interaction between magnons in 
the longitudinal mode. 

We now turn our attention to the difference in the excitation spectra between 
the triangular and square lattice models. This can be easily demonstrated by the 
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Figure 2: The excitation spectrum E(q) of the longitudinal mode together with 
spin- wave excitation spectrum £(q) along (LMTKPQMO) of the BZ with an 
anisotropy A = 1 — 1.5 x 10~ 4 . The gap value at K and Q points is 0.203,25 J for 
the longitudinal mode, and 0.00752s J for the transverse spin-wave modes. 

linear SWT for both models. Due to the frustrations in the triangular lattice, there 
are three sublattice magnetizations with 120° different orientation to one another 
rather than the two sublattices with complete opposite orientation (180°) in square 
lattice. The anisotropy parameter A < 1 is introduced in the y components of the 
Hamiltonian of Eq. (3) rather than in the usual the ordered direction of the z com- 
ponents for the square lattice model, hence the spin- wave spectra of the triangular 
anisotropic model are still gapless at the T point. Furthermore, the anisotropy 
parameter A inside the square root formula of the spin-wave spectrum of Eq. (8) 
is linear rather than the usual quadratic form for the bipartite models such as the 
square lattice model. All these features of the triangular lattice model contribute 
to the slightly softer nature of the longitudinal model than that of the square lat- 
tice model. In particular, the energy gap values for the longitudinal mode is about 
0.203zsJ at Q and at T for the triangular lattice model with the anisotropy 
A = 1 — 1.5 x 10~ 4 as given before, comparing with the bigger energy gap values 
of 0.76zs J at Q and OAAzsJ at T for the square lattice model at the similar value 
of the anisotropy A = 1 + 1.5 x 10 -4 which enters the Hamiltonian in the ordered 
z component [38]. More significantly, due to the fact that antiferromagnetic order 
is noncollinear in triangular lattice, in additional to the above discussed longitu- 
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dinal model corresponding to the ^-component order parameter fluctuations, we 
may have another longitudinal mode corresponding to the perpendicular x-axis 
fluctuations. We leave further discussion to the following section. 

4 Summary and discussion 

In this paper, we have applied our microscopic analysis based on the magnon- 
density waves to study the longitudinal excitations of the z-component order pa- 
rameter fluctuations for a triangular lattice antiferromagnets, with the results simi- 
lar to those of the square lattice model. In particular, the longitudinal spectrum for 
the isotropic model has zero modes at the zone center and the magnetic wavevec- 
tors due to the slow, logarithmic divergence of the structure factor. Also similar 
to the square lattice model, any finite size effect or anisotropy will induce large 
energy gaps at the magnetic wavevectors when compared with the counterparts of 
the transverse spin-wave excitations. Furthermore, we find the longitudinal modes 
of the triangular lattice models are in general softer than that of the square lattice 
model due to the frustrations in the triangular lattice model. Due to the fact that 
the magnetic order in the triangular lattice is noncollinear, there may exist another 
longitudinal mode, corresponding to the order parameter fluctuations in the per- 
pendicular x direction. Instead of the density operator of Eq. (10), we use the 
operator X q = 4 J2i e 4q ' n sf , for the density fluctuations in the x-direction. We 
will report the results for this longitudinal mode elsewhere. 

One remaining question is the intrinsic lifetime of the longitudinal modes for 
triangular lattice model. As the results in Sec. 3 demonstrate that the energy of 
the longitudinal mode is always higher than that of the transverse spin- wave exci- 
tations, we expect the longitudinal mode in the pure triangular lattice model may 
be unstable against decaying into two or more magnons, as emphasized in Ref. 28 
and in a more recent discussion [40]. One way to find the stable longitudinal 
modes is to consider some frustrated models in triangular lattices near a quan- 
tum critical point where the sublattice magnetization is much reduced [41, 42]. A 
more realistic way is to consider some quasi- Id hexagonal systems as mentioned 
in Sec. I. It is therefore particularly interesting to extend our present analysis to the 
hexagonal ABX 3 -type antiferromagnets with both spin quantum number s = 1 
(CsNiCl 3 and RbNiCl 3 ) [23, 24] and s = 5/2 (CsMnI 3 ) [25, 26], where the longi- 
tudinal energy gaps were first observed at the magnetic wavevectors. The spin lat- 
tice structure on the basal planes of these systems are triangular as that discussed 
here. The field theory approach proposed by Afflect based on Haldane's theory of 
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spin- 1 (or integer spin quantum number) chain with five fitting parameters is able 
to explain many features observed, but there are still some disagreements partic- 
ularly for the data away from the minimum energy gap at the antiferromagnetic 
wavevectors [26]. We believe that the longitudinal excitations of the quasi- Id an- 
tiferromagnets is more general than the extension of Haldane's theory of integer 
spin quantum chains and that our microscopic analysis based on the magnon- 
density waves as presented here may provide a more general description of such 
longitudinal modes without any fitting parameter other than those in the model 
Hamiltonian. 
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